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Abstract. - We consider the motion of semiflexible polymers in double-well potentials. We 
calculate shape, energy, and effective diffusion constant of kink excitations, and in particular 
their dependence on the bending rigidity of the semiflexible polymer. For symmetric potentials, 
the kink motion is purely diffusive whereas kink motion becomes directed in the presence of a 
driving force on the polymer. We determine the average velocity of the semiflexible polymer 
based on the kink dynamics. The Kramers escape over the potential barriers proceeds by nucle- 
ation and diffusive motion of kink-antikink pairs, the relaxation to the straight configuration by 
annihilation of kink-antikink pairs. Our results apply to the activated motion of biopolymers 
such as DNA and actin filaments or synthetic polyelectrolytes on structured substrates. 



Introduction. - The Kramers problem [1] of thermally activated escape of an object over a 
potential barrier is one of the central problems of stochastic dynamics. It has been extensively 
studied not only for point particles [2] but also for extended objects such as elastic strings 
which occur in a variety of contexts in condensed matter physics such as dislocation motion 
in crystals [3] , motion of flux lines in type-II superconductors [4] , or charge-density waves [5] . 
Elastic strings activate over potential barriers by nucleation and subsequent separation of 
soliton-antisoliton pairs which are localized kink excitations [6,7]. An analogous problem is 
the activated motion of a flexible polymer over a potential barrier [8] . 

However, the thermally activated escape of a semiflexible polymer, which is a filament 
governed by its bending energy rather than entropic elasticity or tension, remained an open 
question that we want to address in this paper. Semiflexible polymers such as DNA or actin 
filaments have a large bending stiffness and, thus, a large persistence length, L p . On scales 
exceeding L p , the orientational order of the polymer segments decays exponentially, and the 
polymer effectively behaves as a flexible chain with a segment size set by L p . In contrast, on 
length scales which are small compared to L pi the bending energy of the semiflexible polymer 
strongly affects the behaviour of the polymer. The persistence lengths of the most prominent 
biopolymers range from 50nm for DNA [9], to the 10/mi-range for actin [10,11] or even up to the 
mm-range for microtubules [11] and becomes comparable to typical contour lengths of these 
polymers. Whereas the adsorption of such semiflexible polymers onto homogeneous adhesive 
surfaces has been studied previously in [12-14], much less is known about the behaviour of a 
semiflexible polymer adsorbed on a structured surface. 
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In this article we focus on the escape of a semifiexible polymer over a translationally invari- 
ant potential barrier as shown in fig. ^ which can be realized on chemically or lithographically 
structured surfaces. The behaviour of semifiexible biopolymers on such structured substrates 
is of interest, e.g., for electrophoresis applications [15]. Another important class of semifiexible 
polymers are synthetic polyelectrolytes, whose self-assembly and dynamic behaviour on struc- 
tured substrates has only been studied recently [16]. In this article we consider homogeneous 
driving forces across the potential barriers as they can be easily realized on structured sub- 
strates by electric fields for charged polymers as in electrophoresis or by hydrodynamic flow. 
Alternatively, escape over a barrier can be driven by entropic forces arising from asymmetric 
shapes of the potential valleys [17]. 



Fig. 1 - Typical conformation of a semifiexible polymer (thick line) with a kink-antikink pair in a 
double-well potential V which depends on the coordinate z and is independent of the coordinate x. 

Our main results are as follows. As for flexible polymers, the activated dynamics of 
semifiexible polymers is governed by the nucleation of localized kink-like excitations shown 
in fig. ^ We find, however, that the activated dynamics of semifiexible polymers is different 
from that of flexible polymers as kink properties are not governed by entropic elasticity of the 
polymer chain but rather by the bending energy of the semifiexible polymer. This enables 
us to determine the persistence length from kink-properties. Furthermore, we calculate time 
scales for barrier crossing and the mean velocity of the semifiexible polymer for all regimes of 
driving forces: (i) nucleation and purely diffusive motion of single kinks (ii) nucleation and 
driven diffusive motion of single kinks and (iii) for large driving force dynamic equilibrium 
between nucleation and recombination in a kink ensemble. 

Model. We consider the dynamics of a semifiexible polymer in 1+1 dimensions in a 
double-well potential that is translationally invariant in one direction, say the rc-axis as in 
fig-ffl The semifiexible polymer has a bending rigidity k and persistence length L p = 2k /T 
where T is the temperature in energy units. We focus on the regime where the potential 
is sufficiently strong that the semifiexible polymer is oriented along the rr-axis and can be 
parameterized by displacements z(x) perpendicular to the x-axis with —L/2 < x < L/2, 
where L is the projected length of polymer. The Hamiltonian of the semifiexible polymer is 
given by 



V 




z 




(i) 
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i.e., the sum of its bending and potential energy. We consider a piecewise harmonic double- 
well potential V(z) = |Vo(|-2| — a) 2 — Fz that is independent of x and thus translationally 
invariant in the x-direction, where Vo is the depth of the potential and F an external driving 
force density. Below the critical force F c = aVo the potential has two minima at z^ lin = 
±a + F/V a . The Hamiltonian JTJ can be made dimensionless by measuring energies in units of 
a characteristic energy E sc = a 2 n 

1/4^3/4 

, the x-coordinate in units of a characteristic length 
x S c = (^/Vb) 1 / 4 and the z-coordinate in units of a. 

We consider an overdamped dynamics of the semiflexible polymer with an equation of 
motion 

7$* = -^j+C(M) =-Kdtz-V'(z) + a^t) (2) 

where 7 is the damping constant and £ is a Gaussian distributed thermal random force with 
(0 = and (((x, t)C(x',f)) = 2 1 T8{x - x')5{t - if). 

Static kink. At first we construct the static kink for F = 0, which is a localized 
metastable excitation. The static kink z k (x) is defined as the configuration that minimizes 
the energy QJ, i- e -j i s a time-independent solution of (J2J in the absence of thermal noise 
(C = 0) for boundary conditions Zfc(±Z/2) = ±a and d x z k \±L/2 — 0. For F = the potential 
is symmetric and V(z) — V{— z) such that the kink configuration is anti-symmetric with 
Zk(x) = —z k (—x) and centered at x — (i.e. z k (0) = 0). For our piecewise defined potential 
we have to fulfill five matching conditions at x = which connect the two parts x < and 
x > of the kink: z fc (-0) = z k {+0) = 0, d"z k \^ Q = d^z k \ +0 for n = 1,2,3. Both parts 
Zfe(x) + a for x < and z k (x) — a for x > of the static kink are linear combinations of 
the four functions e ± x / w k e ±ix/w k w j iere eight linear expansion coefficients are determined 
from the boundary and matching conditions. The width w k of the kink and the energy E k of 
a single static kink in the thermodynamic limit of large L are given by 

w k = V2x sc = V2(k/Vo) 1/4 and E k = E sc /V2 = a 2 K 1/4 V^ /4 /V2 . (3) 

We expect our results for the kink energy E k ~E sc and width w k ~x sc to hold for all po- 
tentials with a barrier height ~ Vqcl 2 and potential minima separation ~ a independent of 
the particular potential form; only numerical prefactors will differ. We want to point out 
that measurements of the kink width w k and the critical force density F c or the kink energy 
E k are sufficient to determine the bending rigidity n — F c w\/Aa — E k w k /2a 2 and thus the 
persistence length L p = 2k/T if the distance 2a between potential minima is known. 

A static single kink in a polymer of length L is equivalent to one half of a symmetric 
kink-antikink pair configuration with kink-antikink separation d = L in a polymer of length 
2L, as shown in fig.Q] The kink-antikink interaction energy Ei nt (d) = 2(E k (d) — E k (oo)) can 
thus be found by determining the single kink energy in a polymer of length L = d. For large 
separation d/w k 3> 1 we find an exponential decay Ei nt (d) ~ exp(— d/w k ). 

A semiflexible polymer will stay localized to the potential wells even if we set V(z) = for 
\z\ > 2a as long as Vo > Vo, c with Vb, c a 2 ~ (T/L p )(L p /a) 2 ^ 3 according to [14]. This condition 
is equivalent to E k 3> T and thus a small density of thermally induced kink excitations. A 
small kink density in combination with the condition L p ^> a also ensures that the semiflexible 
polymer stays oriented along the cc-axis such that the Hamiltonian JTJ stays valid. The 
condition E k 3> T of a small kink density is equivalent to L p ^> w\j a 2 . For sufficiently strong 
substrate potentials this gives a much wider range of applicability of the Hamiltonian than 
in the absence of a potential where the condition L p > L of weak bending has to be fulfilled 
for a semiflexible polymer to be oriented. 
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Moving kink. - A driving force density F leads to an asymmetry in the potential and 
an effective force on kinks. Moving a kink by — Aa; increases the polymer length in the lower 
potential minimum by Aa: and leads to an energy gain — 2aFAx and thus a constant force 
Tk = —2aF on a kink. As argued above deviations from kink interactions are exponentially 
small for separations d^> Wk- The force Tk leads to kink motion such that we also have to 
consider moving kink solutions. For constant kink velocity v the kink configuration assumes 
a form Zk(x, t) = Zk{x — vt) that solves (0 for ( = 0. Introducing the coordinate y = x — vt 
for the comoving frame, equation (0) reduces to 



ndyZk - vjdyZk + V'(z k ) = 



(4) 



which has to be solved with boundary conditions analogously to the static kink. However, 
in the asymmetric potential the kink is no longer symmetric but centered at yo ^ with 
z k(yo) = where we also have to evaluate the matching conditions. Eq. can be made 
dimensionless by measuring time in units of a characteristic time t sc = 7/V0 and velocities 
in units of v sc = x sc /t sc — k 1 / 4 !/^ 4 /j. For a moving kink both parts Zk(y) — z mm for 
y < and Zk(y) — %rn,in for y > are linear combinations of four functions e KnV where K n 
(n = 1, ...4) are the four roots of the equation kK^ — v^jK n + Vq = that real part of which 
determine the width of the kink w k {v) ~ l/\Re(K n )\. We find K n w k = ±H 1 / 2 (v)±{~H(v)± 
2 3/2 3" 3/ %ff -1 / 2 («)) 1/2 (the first and third sign have t o be id entical) where v = 3 3/4 v/4v sc 
is a dimensionless velocity and H(v) — 3 _1 / 2 ((u 2 + \Jv A — l) 1 / 3 + (v 2 + ^/v 4 — l) -1 / 3 ) an 
increasing, real function with H(v) > H(0) = 1 and H(v) ~ v 2 ^ for v ^> 1. The width of 
the moving kink decreases with velocity as Wk(v) = WkH^ 1 / 2 ^) (for v < 1). In the limit of 
large polymer length L Wk(v) a moving kink solution, fulfilling all boundary and matching 
conditions, must satisfy the force-velocity relation 



F(v) = -F c v 



3 l/4 2 -l/2 g 3/2 ( ^ 

H 3 (v) + 3~ 3 / 2 v 2 



(5) 



For small force densities, we find a linear relation F w — 3 1 / 4 2 1 / 2 F' c u, close to the critical 
force density F c the velocity diverges as — v ~ (1-F/F C )- 3 / 2 , see fig. El 




Fig. 2 - Force density F (in units of F c , solid line) and friction constant r\k (in units of 2>F c a/2 z ^ 2 v s 
3a 2 7/wfe, dashed line) as function of velocity v = 3 3//4 w/4i; sc for a moving kink. 



The result JSJ can also be used to obtain the friction constant r)k of a moving kink by 
equating the friction force vr\k with the driving force Tk = —2aF which gives the relation 
r/fc = 2a\F(v)\/v, sec fig. El % is also related to the energy dissipation rate dE/dt due to kink 
motion which is defined as the product of the friction force — vrjf, and velocity: dE/dt = —v 2 r\k- 
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On the other hand, dE/dt can be calculated directly using the equations of motion J2J) and 
Q in the limit of large L 

— = / dx- — d t z k = -jv / dy{dyz k ) , (6) 
dt Sz k J_ OQ 

and we read off a kink friction constant r)k = 7 f-^ dy (d y z k ) 2 . Integration in the limit of 
small driving forces gives n k ~ 3ja 2 /2w k and equating the friction force with the driving force 
—2aF — vr/k gives a linear relation v = —4Fw k /3ja which agrees to leading order with our 
above result (J5J , see also fig. |21 

Thermal noise and kink motion. For a more detailed analysis of the effect of noise 
on the kink motion we consider noise-induced perturbations of shape and velocity of a kink 
moving with constant velocity v. For a time-dependent kink center at Xk(t) the comoving 
frame coordinate is given by y = x — x k (t) . Adding shape perturbations to the corresponding 
kink solution z k (y) of (J2J, we arrive at the decomposition 

oo 

z(x,t) = z k (x - x k (t)) + ^X p (t)(/) p (x - x k (t),t). (7) 

P =i 

4> p are normal modes of the kinked polymer which we will determine below and X p (t) are 
expansion coefficients; the zero mode of kink translation is explicitly taken into account by 
positioning the kink center at x k (t). Substituting JJJ) into the equation of motion J3J) and 
expanding about the kink, we obtain 

(oo \ oo 

d v z + ^ X p dy(j>p J + j^Xp^p = C (8) 
P =i / p=i 

if the normal modes 4> p (y,i) = / P (y)e~" pt fulfill the eigenvalue equation 

Kd^fp - ivdyfp + V"{z k (y))f p = ujpjfp (9) 

where V"(z) = Vq(1 — 2a6(z)). @ has to be solved with boundary conditions f p (—L/2) = 
f p (L/2) = and f' p {-L/2) = f p (L/2) = where we consider the limit L/2 > x k (t) and 
neglect the shift of boundaries in the comoving frame. Eq. 10 has a set of eigenvalues ui p with 
orthonormal eigenfunctions f p (y) (with respect to the scalar product (f\g) = J dyf(y)g{y)). 
The translation mode /o = dyZ k (y)/C of the kink has zero eigenvalue ujq — 0. C is a normal- 
ization constant determined by C 2 — {dyZ k \dyZ k ). Multiplying eq. ||SJ| with the translation 
mode fo and integrating yields an equation of motion for the kink 

-l 

(10) 

where Cfe(^) = ^i^l)^ 1 J dyfo(y)£(y + x k (t),t) is an effective Gaussian thermal noise for 
the kink with correlations (Cfe (*)Cfe (*')> = 5 ( l - *')( 2 r/C 2 7) (where we used (/olio) = !)■ 
The sum in (f 101) represents terms from kink-phonon scattering neglecting of which leads to 
an overdamped Langevin equation x k {t) = v + Ck(t) describing Brownian motion with drift. 
From the noise correlations we can read off the corresponding diffusion constant of the kink 
as D k = T/C 2 "f. Note that the corresponding kink friction constant r\ k = T/D k is identical 
to our above result JSJ obtained from complementary energetic considerations. 



x k =v + Cfe 



l + c- 1 yZXpe-^iMdyfp) 



p=i 
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Fig. 3 - An ensemble of well-separated kinks and antikinks which move with velocity v and — v, 
respectively. 



If kink-phonon scattering is neglected, the kink is performing a Brownian motion with 
drift. The polymer crosses the potential barrier by moving a kink over the entire length L of 
the polymer. Thus, the average crossing time is t cr ~ L/v for the case of directed diffusion 
with v > under the influence of a driving force density F. For F = and v — the kink 
performs an unbiased random walk with (x 2 ) ~ D^t from which we estimate the average 
crossing time as t cr ~ L 2 rjk/T ~ L 2 ja 2 /Twk- For F = and at low temperatures t cr gives 
the relaxation time from a kinked state as in fig. ^ to a kinkless state. The diffusive part of 
the kink motion can be neglected for forces F 3> 2T/La. 

Thermally activated barrier crossing. For sufficiently large F, thermally activated 
barrier crossing proceeds by the nucleation and subsequent separation of a kink-antikink pair, 
see fig. |3 Each passing kink or antikink increases the polymer position by Az = 2a. For an 
ensemble of pL kinks and pL antikinks with kink density />< 1 /wk as in fig. the fraction 
of moving polymer segments is given by 2pLutk/L = 2pwk- These polymer segments move 
with velocity 2av/wk in the z-direction which leads to the average velocity v z = (dtz) — Aavp. 
The kink density p is determined by the dynamical equilibrium of kink nucleation with rate 
j (per length) that we will calculate below, see l|llf) . and kink-antikink recombination with 
rate 2p 2 v [6]. Equating both rates gives a steady-state density p 2 = j/2v and thus an average 
polymer velocity v z = 2a(2vj) 1 / 2 . 

In order to find the nucleation rate j, we use Kramers theory. In the following, we only 
give the main results of this calculation, details will be described elsewhere. As for flexible 
strings [6] , the dynamics of the nucleation is governed by the critical nucleus representing the 
saddle point in the multi-dimensional energy landscape. The critical nucleus configuration 
z n (x) is the analogon of the static kink-antikink pair for F > and fulfills the same saddle- 
point equation 8TL/8z — 0, see (J2J- For the critical nucleus we obtain an excess energy 
AE n sw 2Ek(l — F/F c ) 2 that enters the nucleation current j ~ exp(—AE n /T). The prefactor 
depends on the corresponding attempt frequencies and, thus, the spectra Lo n ,p and lj s>p of 
fluctuations around the the critical nucleus and the straight polymer, respectively. For the 
straight polymer we find a spectrum of stable phononic modes LJ S fl — V0/7 and oj s ^ p w 
V0/7 + k((tt/2 + pir)/L) 4 /j (p > 1). For the critical nucleus, the spectrum consists of an 
unstable mode w n ,o < representing the collective coordinate along which the nucleation 
proceeds, a zero mode Lo n ^\ — corresponding to the translation of the nucleus, and a sequence 
of stable phononic modes uj n ,2 = V0/7 and u> ntP ~ Voh + n{{b+ P Ti)/L) i / 1 (p > 3), where 
& is a numerical constant. Using Kramers theory in the regime F > T/2awk [18], we finally 
obtain the nucleation rate 

j = {2^ 2 1 1 l 2 GT- 1 l 2 Q n exp (-AE n /T) (11) 

where Q 2 n = \uo nfi \uj sfi uj s s]\ p>1 uo s ^/uj n ^ w |1- 2 4 / 3 (l - F/F c y 8 / 3 \(V /"/) 3 contains all at- 
tempt frequencies, and G = L^ 1 J dx n [J dx (d x z n (x)) 2 ] 1 ^ 2 w a(l — F / F c ) / ^/vJk is the Jacobian 
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for the change of coordinates from the amplitude of the translational mode d x z n to the nucleus 
position x n . 

For small driving force densities F <C 2pT/a, the kink motion is diffusive, and the above 
approach breaks down as kink-antikink pairs cannot separate but tend to recombine. For 
Ek 3> T, the system reaches thermodynamic equilibrium with a low kink density p eq ~ 
exp (— Ek/T) given by the Boltzmann distribution and with v z = <lavp e q. For intermediate 
driving forces IpTja <C F < F/2awk, the critical nucleus is in quasi-equilibrium [18], and 
we find again j ~ exp(— AE n /T) as in the high-force expression l|lll) but with a different 
parameter dependence of the prefactor. 

Conclusion. - In conclusion, we described the activated dynamics of semiflexible poly- 
mers which is governed by kink excitations. The static kink has the energy E^ and the width 
uik as given by © . Both kink properties are governed by the bending rigidity of the semiflex- 
ible polymer. In the presence of a driving force density F, there is a force Tk acting on the 
kink that leads to moving kink solutions which satisfy the force- velocity relation (JjJJ. In the 
absence of kink-phonon scattering the kink performs Brownian motion with drift for which 
we have calculated the friction constant r\k and the diffusion constant Dk- This leads to es- 
timates for the crossing times t cr ~ L/v for large forces F ^> 2T/La and t cr ~ L 2 rjk/T for 
small forces F <C 2T/La. For large forces, the nucleation of kinks proceeds by activation over 
the saddle point which represents the critical nucleus. Application of Kramers theory leads 
to the nucleation rates (fill) which determine the average velocity (d t z) of the polymer. Our 
results are not only relevant to the dynamics of semiflexible polymers but can be extended to 
kink excitations in fluid membranes [19]. 
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